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TIDAL DECAY AND CIRCULARIZATION OF THE ORBITS 
OF SHORT-PERIOD PLANETS 

Adrian Rodriguez 1 and Sylvio Ferraz-Mello 1 



Abstract. We analyze the long-term tidal evolution of a single-planet 



system through the use of numerical simulations and averaged equa- 
tions giving the variations of semi-major axis and eccentricity of the 
p£ ■ relative orbit. For different types of planets, we compute the variations 

due to the planetary and stellar tides. Then, we calculate the critical 
value of the eccentricity for which the stellar tide becomes dominant 
l—i ' over the planetary tide. The timescales for orbital decay and circular- 

' ization are also discussed and compared. 

, 1 Introduction 

OH It is well-known that tidal friction produces variations in the orbital elements of 
a close-in companion. In the case of an interacting pair with two extended (and 
tidally deformed) bodies, the tidal evolution depends on the rotational state of the 
bodies. We can identify two important examples: the satellite and the exoplanet 
£Q cases. In the first case, the rotation of the central body (a parent planet) is 
generally much faster than the mean orbital motion of the satellite orbiting the 
| primary. In the exoplanet case, the rotation of the central body (a parent star) 
is generally much slower than the mean orbital motion of the planet. Changes 
in orbital elements are also accompanied by variations in the rotation of each 
deformed body. The reader is referred to Ferraz-Mello et al. (2008) for further 



o 



H ' details. 



2 Mean orbital variation 

Consider a two-body system formed by a single star and a short-period companion 
planet. We suppose that both bodies are able to be tidally deformed due to the 
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mutual interaction. Tides on each body provokes variations in the elements of 
the relative orbit and their rotations. The equations that govern the mean orbital 
changes in the astrocentric orbital elements of an exoplanet are, following Ferraz- 
Mello et al. (2008) (correcting a misprint in equation (93) and neglecting the 
inclinations and the terms proportional to O/n): 

(n) = ^%^4^(l + 23e 2 + 7£>e 2 ) (2.1) 

(e) = _ Wnek d *m p Rl^ I i + lj I ,2.2) 




2m,o 5 



where n, a and e are the mean orbital motion, semi-major axis and eccentricity 
while m* ,R* , kd* and m p ,R p , kd P are the masses, radii and dynamical Love numbers 
of star and planet, respectively. We stress that the above equations are valid only 
in the case of star-exoplanet interacting pair. The main equations for the planet- 
satellite case are different (Ferraz-Mello et al. 2008). 
The parameter D is defined as 




D = ^p_^i — ir ( 2 - 3 ) 



where e^* and e' 2p are lag angles associated to the tidal waves whose frequencies 
are 217* — In (on the star) and 2Q p — n (on the exoplanet), respectively. Here, Q is 
the angular velocity of rotation of the tidally deformed body. The lag angles come 
from the delay in their response to the tide raising potentials (see Ferraz-Mello et 
al., 2008 for details). Several hypotheses were done to obtain the above equations. 
The planetary rotation was assumed in a quasi- synchronous state (f2 p ~ n) and 
the star rotation verifies f2* <C n. For the star, a linear model is assumed in the 
relationship between lag angles and frequencies (Darwin 1880; Mignard 1979). In 
this case, each lag angle is proportional to the frequency of the corresponding tidal 
wave, and the coefficient of proportionality At* (time delay), is the same for all 
frequencies. No model needs to be assumed for the planetary tide lags, except that 
equal frequencies are assumed to span equal lags. The motion is supposed planar, 
i.e the reference and orbital planes coincide (zero obliquities). The equations are 
valid only up to second order in the eccentricity. 

In order to simplify the notation, we define the following parameters: 



9fcd*mp^ 5 „ _ 9k dp m 

4 Q* m* 2 Q p m 

e^- 1 and Q p = \e' 2p \~ l 
kd and Q are poorly known quantities for both stars and planets). Hence 



- 7 7f— ; /- -7T 1 — ( 2 - 4 ) 

where Q are quality factors defined by Q* = |eg„| -1 and Q p = le^l -1 (in general 
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Introducing (|2.4|) into equations (|2.ip and (|2 . 2[) and using the third Kepler law to 
relate a to h, (a) = — ^(n), we have 

(d) = -|na~ 4 [(2 + 46e 2 )s + 7e 2 p] (2.6) 
(e) = -inea~ 5 (18J + 7p). (2.7) 

o 

Note that, for each equation, the terms proportional to s andp are the contribution 
to the total variation of the tides raised on the star and planet, respectively. 

3 Timescales of tidal evolution 



It is clear from equations (|2.6[) and (|2.7p that the effect on the orbital elements 
due to the tidal interaction is to reduce the size and eccentricity of the relative 
orbit. In this section, we discuss the decreasing timescales of semi- major axis and 
eccentricity. 

3.1 Semi-major axis 

The decreasing timescale of the semi-major axis can be defined as r a = a/|d|. 
Using equation (|2.6p . 

r a = ^n~V[(2 + 46e 2 )J + 7e 2 p}-\ (3.1) 

If we need to know the timescale only due to planetary tides, it is enough to put 
s = (or Q^ 1 = 0) in the above equation to obtain 

Note that limr^ = oo as e — > 0, indicating that, when the only tide is the planetary, 
the semi-major axis stops decreasing after circularization. The contribution of the 
stellar tide follows a similar analysis putting p = (or Q" 1 = 0). The result is 

a ~ 2(2 + 46e 2 )J- 1 j 

e -> 0. This 

the semi-major axis continues to decrease due to the stellar tide 



Note that limr* = 3n ,. a < oo as e — > 0. This shows that after circularization 

' 4s 



3.2 Eccentricity 

The timescale of orbital circularization can be defined as r e = e/|e|. Using (12. 7p 

Te = m^7-p (3 - 4) 
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As in the case of the semi-major axis, we can compute the individual timescale 
due to each tide. For the planetary and stellar tides we have, respectively 



3n 1 a 5 „, 3n 1 a 5 



7p ' e 185 
Note that rf and r* are independent of e. Moreover, we see that 

Tj_ = 18 £ 
T* 7 p" 



(3.5) 



(3.6) 



In order to have an idea of the numerical values of the above quantities, we give 
a simple example. Consider a planet-star system with masses and radii equal to 
those of the Sun and Jupiter. We have 

|^ 50 |dpQ* (37) 

In the above equation, typical values are 2.5 x 10 4 < Q p /kd p < 2.5 x 10 5 (Laskar 
& Correia 2004) and 2 x 10 7 < Q*/fc d * < 1.5 x 10 9 (Carone & Patzold 2007). 
Hence > 1 and § > 1. Then from t£ < r*, meaning that the tidal 

evolution of the eccentricity is mainly determined by planetary tide. In the case 
of a Super-Earth (m p ~ 5m®), we have fc fc d j^ ^> 1 because the quality factor 
of a rocky planet is several orders of magnitude smaller than the corresponding 
value for a gaseous planet. In this case, the stellar tide is completely negligible 
compared with the planetary tide in producing orbital circularization. 

4 Orbital decay for different type of planets 

4.1 Critical eccentricity for orbital decay 

As we have done in the case of timescales for orbital circularization, we may 
compare the timescales for orbital decay due to planetary and stellar tides. For 
that sake, using equations (|3.2[) and (|3.3|) we obtain 

t* 7e 2 p 

Note that the above timescales ratio depends on the eccentricity, which is a time- 
dependent quantity. We can compute a critical value of e for which equation (14. ip 
is equal to unity. It is 



7p/s - 46 



(4.2) 



If e > e c then < t* and the planetary tide dominates. On the contray, if e < e c 
then tP > t* and the stellar tide becomes more effective to produce orbital decay. 
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Fig. 1. A plot of equation 114.21 1 with p/s as independent variable. See that lime c = as 
p/s — > oo. 



Figure [T] is a plot of e c as a function of p/s. The minimum value of p/s for 
which equation (|4.2| has a possible solution (e c < 1) is 48/7. Note that, for 
high values of p/s, it is sufficient a small value of eccentricity for that the orbital 
decay is dominated by the planetary tide (tJ < r*). On the other hand, for very 
small values of p/s, an eccentricity close to unity is necessary to have comparable 
evolutions due to planetary and stellar tides (remember that the model is valid up 
to second order in e). 



4.2 Hot Jupiters and hot Neptunes 

In this section we analyze the tidal evolution considering two different types of 
close-in planets, depending on the planetary mass. In a general way we call hot 
Jupiter any close-in planet for which m p ~ mj, and hot Neptune planets for which 
m p ~ mjv, being mj and mjsi the mass of Jupiter and Neptune, respectively. From 
equations lj3.1|) and (|3.4|) we know that 



Tg = 1 18s + 7j3 

r e 2 [(2 + 46e 2 )J + 7e 2 P ] ' 1 ' j 

It is clear from the above equation that, for small values of e, r a ^> r e . This indi- 
cates that the effect of tides circularize the orbit before producing any appreciable 
decrease in semi-major axis. Note that this fact minimizes the contribution due to 
the planetary tide because after circularization, only the stellar tide contributes to 
produce the orbital decay (see equation (|2.6[) V Equation (|2.4p shows that the stel- 
lar tide contribution is proportional to the planet mass m p , where the coefficient 
of proportionality depends on stellar paramaters only. Hence, massive planets are 
more effective to have orbital decay due to the stellar tide, indicating that hot 
Jupiters must decay faster than hot Neptunes. 



6 



Title : will be set by the publisher 



Planet 


m p (mj) 


Qp 


a (AU) 


eo 


p/s 


e c 


Jupiter-like 


1.0 


1 x 10 b 


0.04 


0.1 


125 


0.05 


Neptune- like 


0.1 


1 x 10 4 


0.04 


0.1 


2.40xlO a 


0.01 



Table 1. The data for numerical simulations. Note that the value of ao corresponds to 
an initial orbital period of Po = 2.92 days. The radii values are computed imposing that 
~p p = P/jv- For the star m* = rag, R, = Rq and Q» = 1 x 10 6 . 



5 Numerical simulations 

In this section, we give some examples corresponding to the results obtained in 
the preceding sections. We adopt a linear model (Mignard 1979) to relate lags 
and frequencies of tidal waves for the star and the planet. It is important to 
mention that the Mignard's model gives the tidal force in a closed form, without 
any expansion in eccentricity. In that model the lag angles are given by 6j = z/,At, 
where is a tidal frequency and At is a measure of the delay in the response of 
the tidal deformation, which is assumed equal for all tidal waves. It is easy to 
make the passage from Q to At provided the tidal frequency is fixed and taking 
into account that Q = |e| . We then numerically compute the evolution of a 
single-planet system composed by a Sun-like star and a short-period companion 
planet. The selected input parameters are given in table [1] for two different types 
of planets. In each case, the values of p/s and e c are computed through equations 
(12. 5|) and ()4.2|1 . For all examples showed in this section we have included a scaling 
factor of 10 4 multiplying the tidal parameters of both bodies in order to accelerate 
the computation. 

Figure [2] shows the long-term tidal evolution of semi-major axis and eccentricity 
in the case of a hot Jupiter around a Sun-like star. One may note that the timescale 
for orbital circularization r e is much smaller than the timescale for orbital decay 
T a , as discussed in section 1431 The figure also shows the comparison of the cases 
with and without stellar tide. It is well appreciable how the inclusion of the stellar 
tide changes the evolution of the semi-major axis, compared with the case where 
only the planetary tide is considered. Figure [5Jd is a zoom of and it shows the 
time from which the stellar tide becomes dominant. Comparing to figure [2]:, we 
note that this time corresponds to the moment in which e ~ e c = 0.05. When 
e > 0.05, the planetary tide dominates. However, since rf <C t%, the orbit reaches 
very quickly the critical eccentricity and after that the stellar tide starts to be 
dominant. We note that the eccentricity tidal evolution is almost independent on 
the inclusion of the stellar tide. This happens because rf <C t* and therefore the 
contribution of the stellar tide to the orbital circularization is almost negligible. 

Figure [3] shows the long-term tidal evolution for a Neptune-like planet inter- 
acting with the same star under the same initial conditions (m p ~ 2mjv). As in 
the case of a Jupiter-like planet, the semi-major axis evolution is different when 
the stellar tide is added to the planetary tide. However, in this case, the difference 
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0.05 0.1 0.15 0.2 0.25 0.3 
time [Gyr] 

Fig. 2. (a) Long-term tidal evolution of semi-major axis for a Sun-Jupiter tidal interacting 
system. The cases with and without the stellar tide are showed in order to illustrate the difference 
between them, (b) Zoom of figure (a) showing the moment when the planetary tide becomes 
smaller than the stellar tide because the orbit reaches the critical eccentricity, (c) Eccentricity 
tidal evolution. The curves corresponding to the planetary and total tide are superimposed (see 
text for discussion). 

is not so apparent as in the case of a Jupiter-like planet, at least for the same 
timescale (compare figure with [2k). Note that r a jv ~ 10r a j. On the other 
hand e C M — \ e cj, meaning that the planetary tide dominates almost until total 
circularization. After that point, the stellar tide controls the orbital decay. More- 
over, since -C mj, the influence of the stellar tide in the case of Neptune- like 
planets is not so effective in producing strong damping in semi-major axis as in the 
case of a Jupiter-like planet. This was already anticipated in section l4~2l where we 
have seen that s oc m p . The eccentricity evolution does not show any appreciable 
change when the stellar tide is included, as we see in figure [3)d. This is because 
the relation rf <C r* is more closely satisfied for the case of Neptune-like planets 
(see section I3.2I) . Note also that r e jv < r e j . 

6 Conclusion 

The main conclusion of this communication is that, in terms of timescales, the 
stellar tide plays an important role to determine the orbital decay of hot Jupiter 
planets. In the case of hot Neptune- like planets the contribution of the stellar tide 
is smaller but still important. The timescale for orbital decay is not sensitive to the 
value of Q p , since the long-term tidal evolution of semi-major axis is controlled 
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Fig. 3. (a) Neptune-like planet semi-major axis tidal evolution. See that the small planetary 
mass obliges to a much slower semi-major axis decay compared with the Jupiter case. It is 
important to note that T a is higher than the age of the universe, (b) Zoom of figure (a).(c) 
Eccentricity tidal evolution. Note the moment for which e = e c = 0.01. 



by the stellar tide. The orbital circularization is well-determined only by the 
planetary tide and is Q p dependent. The computation of e c enable us to known 
the moment for which the stellar tide becomes dominant. 
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